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Введение. Исследованы резонансные частоты продольных 
колебаний системы, состоящей из цепочки разных по геомет- 
рическим и физическим параметрам элементов: пьезокерами- 
ческого, упругого и акустического. Проведено сравнение ре- 
зультатов в пакетах СОМ$ОГ и АСЕГАМ. Оценена зависи- 
мость значения первой собственной частоты от геометриче- 
ских параметров. Исследовано влияние динамической вязко- 
сти на АЧХ продольных колебаний. 

Материалы и методы. Для исследования выбрана система 
упругих и акустических элементов, позволяющая описывать 
работу ультразвукового режущего устройства. Возбудителем 
колебаний выступает пьезокерамический преобразователь, 
совершающий колебания по толщине. Концентратор колеба- 
ний и стержневой элемент принимаются в качестве упругих 
элементов и выполнены из нержавеющей стали. Имитатором 
режущего элемента выступает акустическая жидкость. Про- 
веден модальный и гармонический анализ сложной системы, 
состоящей из разных по физическим свойствам элементов. 
Результаты исследования. Построены осесимметричная и 
трехмерная конечноэлементные модели исследуемой систе- 
мы. Для концентратора колебаний предложены различные 
виды кривизны и толщины звена с изменяемой формой по- 
верхности. Получены первые собственные частоты продоль- 
ных колебаний стержневого элемента, контактирующего с 
акустической жидкостью. Установлено хорошее согласие с 
результатами работы конечноэлементного пакета АСЕГАМ. 
Получены амплитудно-частотные характеристики колебаний 
концевой части вблизи первой резонансной частоты. Отмече- 
но, что вязкость акустической среды мало влияет на амплиту- 
ду колебаний упругого стержня и совсем не влияет на резо- 
нансную частоту. 

Обсуждение и заключения. Гармонический и модальный ана- 
лиз показал, что высокочастотные продольные колебания 
стержневого элемента значительно зависят от тангенса угла 
потерь упругих элементов и слабо зависят от вязкости контак- 
тирующей акустической среды. Полученные результаты могут 
представлять интерес при конструировании ультразвуковых 
режущих медицинских приборов. 
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Введение. Ультразвуковые скальпели позволяют проводить рассечение и коагуляцию тканей с целью оста- 
новки кровотечения [1-4]. Конструктивно такие приборы состоят из цепочки упругих и пьезокерамических (или маг- 
нитострикционных) элементов, собранных в единую систему, создающую максимальные продольные колебания ре- 
жущему элементу на частотах от 20 до 50 КГц. Наиболее эффективным будет режим, при котором режущий элемент 
совершает продольные колебания на резонансной частоте. Под воздействием кавитации клеточной структуры приле- 
гающих к лезвию слоев ткани сухое трение переходит в полусухое или даже жидкостное. Это приводит к существен- 
ному уменьшению как нормального, так и тангенциального усилия резания, но максимальная амплитуда колебаний и 
частота резонанса изменяются. Поэтому указанная в инструкциях к ультразвуковым скальпелям основная частота ко- 
лебаний без учета нагрузки на режущий элемент вряд ли будет резонансной. С другой стороны, в цепочке упругих 
элементов содержится концентратор колебаний, позволяющий увеличить амплитуду продольных колебаний рассека- 
ющего элемента, форма и размеры которого влияют на амплитудно-частотные характеристики скальпеля. Еще одним 
важным элементом является пьезокерамический преобразователь, который может быть многосекционным, а может 
иметь и неоднородную предварительную поляризацию [5—7]. 

Для моделирования работы ультразвукового режущего прибора была рассмотрена конечноэлементная модель 
составной системы пьезокерамических и упругих элементов, один из которых имеет превалирующий над остальными 
продольный размер, с целью изучения влияния параметров физического и геометрического характера на его продоль- 
ные колебания. Такая система близка к конструкции медицинского ультразвукового скальпеля и может быть исполь- 
зована при моделировании его работы. 

Среди факторов, которые влияют на собственные частоты исследуемых колебаний, можно выделить три ос- 
новные группы: геометрические параметры, материальные свойства, неоднородная поляризация пьезоэлемента. Кро- 
ме этого, при расчете рабочих характеристик учитывалась среда, с которой контактирует рассекающий элемент. С 
этой целью была использована модель вязкой акустической жидкости. Для анализа поведения скальпеля в различных 
типах тканей параметрами этой жидкости (например, вязкостью и плотностью) можно варьировать. Различные подхо- 
ды к моделированию свойств живых тканей, в том числе и акустических жидкостей, представлены в [8, 9]. 

В данной работе построена параметризованная геометрия модели, проведена серия численных экспериментов 
и определены изменения собственных частот в зависимости от размеров отдельных элементов системы. Кроме того, 
были проанализированы амплитудно-частотные характеристики системы на частотах, близких к резонансу в зависи- 
мости от физических свойств акустической среды жидкости. 


Материалы и методы. Для исследования была взята система, состоящая из цепочки упругих и пьезокерами- 
ческих элементов, показанных на рис. 1. 








Рис. 1. Параметризованная геометрия модели скальпеля в осесимметричной и в трехмерной постановке в пакете СОМ$ОГ 
Е р. 1. Рагатетте4 зсаре[ то4е! веотету т амзуттейчс апа йтгее-@тепяопа! ет т СОМ5ОГ расказе 


Здесь участки АВ, ВО, ОС — пьезокерамический элемент, концентратор колебаний и стержневой элемент со- 
ответственно. Все элементы имеют осесимметричное строение. Пьезоэлемент представляет собой цилиндрический 
пьезокерамический преобразователь с толщиной В и предварительной поляризацией, направленной по толщине. Элек- 
троды нанесены на лицевые поверхности пьезокерамического цилиндра, и на них подается разность потенциалов, из- 
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меняющаяся по гармоническому закону во времени. Материал концентратора колебаний и стержневого элемента — 
нержавеющая сталь. Контактное взаимодействие между преобразователем и концентратором колебаний и между кон- 
центратором и стержнем считается безотрывным. Сечение А преобразователя закреплено; на боковой поверхности 
имеет место гладкий жесткий контакт. 

Последнее условие выполнено и для концентратора напряжений на участке ВС длиной а. Остальная боковая 
поверхность длиной Ъ свободна от механических напряжений. Для стержня принималось условие гладкого жесткого 
контакта на участке РЕ длиной с. На участке ЕЁ длиной 4 принимались нулевые напряжения. На участке ЕС длиной е 
и на его концевом сечении С имеет место сопряжение с акустической средой КОММ. При этом считалось, что на ли- 
цевой поверхности, проходящей через точку К, звуковое давление равно нулю, на боковой поверхности КТ, имеет ме- 
сто условие импедансного типа, а лицевая поверхность ММ — непроницаемая плоскость. 

Каждая из отмеченных выше сред описывается линейными моделями, которые для упрощения представим в 
векторном виде. 

Пьезокерамическая среда: 


У.в=ро” и; 6=С:=-е'.Е; ==1/2(Уш+Уи'); 
У.0=0; Ю=е:=+э.Е; Е =-УФ; 


в которой 6,50, Е, и, ф,С,е,э,р,« — тензор напряжений, тензор деформаций, вектор электрического смещения, век- 


(0 


тор электрического поля, вектор перемещения, электрический потенциал, тензор упругих постоянных, тензор пьезо- 
электрических констант, тензор диэлектрических проницаемостей, плотность материала керамики и круговая частота 
соответственно. 

Упругая среда: 


У: 6 =р, 0? и; б=Лх НЕЕ, #; &=1/2(Уш+Уи'); (2) 
в которой в,#, и, Лк, их,Рх — тензор напряжений, тензор деформаций, вектор перемещения, единичный тензор, ко- 


эффициенты Ламе и плотность материала концентратора колебаний и стержня соответственно (А = 1,2). Затухание в 
этих материалах моделируется с помощью комплексных коэффициентов Ламе, и тангенс угла потерь равен 0,0001. 
Акустическая среда: 





Р 0.5 
1 ор [0 | > + с? 
у. Ур Ч 0; = |; в=е ив ь ре = : (3) 
Ре Ре Сс рс Се 


В которой р,р,с,иа.иь —_ акустическое давление, плотность жидкости, скорость звука в ней, динамическая вязкость и 


второй коэффициент вязкости среды соответственно. 

Граничные условия и условия сопряжения были подробно описаны выше, они очевидны, поэтому математи- 
ческие формулы для них выписываться не будут. Добавим только, что, если будут рассматриваться неоднородно по- 
ляризованные керамические таблетки, то упругие, пьезоэлектрические и диэлектрические характеристики таких сред 
будут представлять функции координат области, занятой керамикой. 

Для описания геометрии всех элементов системы были разработаны три конечноэлементные модели. Две осе- 
симметричные модели были построены в пакетах АСЕГАМ и СОМЗОГ. Третья модель построена в пакете СОМ$ЗОГ 
для случая ЗО. 

В основе каждой модели лежит параметризованная геометрия, представленная на рис. 1, где также указаны 
независимые параметры, отвечающие за геометрию модели. В левой части системы расположен пьезокерамический 
цилиндр, выделенный более темным цветом. Следующий элемент — стальной концентратор колебаний, состоящий из 
двух геометрических компонентов: цилиндра и конусообразного переходника. В реальных конструкциях хирургиче- 
ских инструментов толщины преобразователя и радиусы концевых сечений концентратора связаны с длинами распро- 
страняющихся волн. В предлагаемом варианте мы отказались от дополнительных ограничений на геометрические 
параметры, чтобы иметь возможность получить более общие результаты. При этом в модель заложена возможность 
изменять функциональным образом диаметр конусообразной части. С этой целью образующая конуса в глобальных 
координатах в плоскости У, описывается квадратичным уравнением: 


—_ 2 __2 
У(2) = А? +КЕ+К,; бы” ВЫ в. 


2-20 





К, =7-А 20; 2, =йс+а; 2 =2, +В, 
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которое позволяет проводить численные эксперименты для вогнутых и выпуклых профилей направляющей, изменяя 
лишь один параметр К, как это показано на рис. 2. 





Рис. 2. Различные формы переходника между пьезоэлементом и стержнем (форма определяется параметром А) 


Ев. 2. Татоиз югтз ор ше адартег Белееп ше рлего@есичс сеЙ ап4 Ше го4 (Ше зйаре 15 4йегттте4 Бу Ше рагатёег К) 


Следующий в цепочке упругих элементов — стержень, имитирующий тело рассекателя. Стержень разделен 
на три условные части, на каждой из которых задаются свои граничные условия, как говорилось выше. Акустическая 
часть также представлена на рис. 1 затемненным цветом и располагается справа. Для удобства в этой среде был сде- 
лан разрез, позволяющий оценить ее влияние на продольные колебания стержневой конструкции. Полный список не- 
зависимых геометрических параметров приводится в табл. 1. 












































Таблица 1 
Та 1 
Изменяемые параметры модели 
Мо4е срапетз рагатейегу 
Обозначение | Начальное значение Описание 
7. 7,5 мм Радиус керамического диска 
й 5 мм Толщина керамического диска 
а 5 мм Толщина слоя металла 
Ь 10 мм Толщина переходника 
Г, 2,5 мм Радиус стержня 
К —80 Безразмерный параметр кривизны образующей конуса переходника 
С 0,02 Длина начального сегмента стержня 
а 0,2 Длина среднего сегмента стержня 
е 0,03 Длина сегмента стержня, погруженного в жидкость 
Ц Ра-5 Вязкость жидкости 
р 1500 кг/м3 Плотность жидкости 
с. 1900 м/с Скорость звука в жидкости 














Физические параметры жидкости, симулирующей поведение среды, с которой происходит контактное взаи- 
модействие, находятся в нижней части этой таблицы. Это скорость распространения звука, плотность и динамическая 
вязкость (второй коэффициент вязкости принимался равным нулю). Более подробно математические модели электро- 
упругих и упругих тел (1)-(2), нагруженных на акустическую среду (3), используемые в пакете АСЕГАМ, представле- 
ныв [10]. 

Результаты исследования. Для осесимметричной и трехмерной моделей была проведена оценка качества 
конечноэлементных сеток. Задача на собственные значения решалась на сетках с различными размерами. Измельче- 
ние сетки продолжалось до тех пор, пока дальнейшее уменьшение размеров не начинало влиять на результат по за- 
данной погрешности. В нашем случае результатом исследования была выбрана первая собственная частота резонанса 


продольных колебаний |. Продольные колебания отбирались посредством перебора собственных форм колебаний. 


В табл. 2 приведены результаты расчетов, по итогам которых определены оптимальные размеры сеток, позволяющие 
получить достоверные результаты с минимальным числом элементов. 
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Таблица 2 
ТаЫе 2 
Результаты подбора сеток 


Спа з@есНоп гезий5 




















2 ЕВ) 
Сша) | Л 9 | сюбыы | Л 0 | сиболь 
[6,8 -104; 0,034] 4900,8 827 5050 3300 
[1,02:10; 0,0228] 4900,3 1059 5011,4 6316 
[1,02:10; 0,018] 4900,1 1095 4934,7 22960 
[4,25:10-5; 0,0126] 4900,1 1559 4932,1 178816 























Представленные результаты позволяют сделать вывод, что для достижения заданной точности осесимметрич- 
ная модель позволяет использовать значительно меньшее число степеней свободы. Так, с погрешностью б = 0,6% в 
осесимметричной модели можно взять в 115 раз меньше конечных элементов, что важно при решении задач оптими- 
зации по параметрам. Однако в дальнейшем возникнет необходимость использовать трехмерную модель, когда будет 
учитываться неосесимметричный вид наконечника стержня и неосесимметричная форма внешних воздействий. 

Первый этап численного эксперимента заключался в валидации моделей путем сравнения различных реализа- 
ций при решении задачи на собственные значения. Была отработана методика определения первых трех рабочих ча- 
стот устройства, которые соответствуют продольным колебаниям стержня (в табл. 2 указано значение только первой 
из НИХ). 

Следующий шаг численных экспериментов заключался в сравнении результатов расчета собственных частот 
продольных колебаний, полученных разными моделями и разными конечноэлементными комплексами. 
АСЕГАМ была 


табл. 3 представлены результаты вычисления первых двух собственных частот. 


В частности, в пакете построена осесимметричная модель этой системы. В 


Таблица 3 
ТаЫе 3 
Результаты решения задачи на собственные значения на оптимальных сетках в различных пакетах 


Везийу оГ5о[уте евепуаше ргоет оп орнта[ эта; т уапои5 расказез 














Модель р (Ги) . (Гц) 
АСЕГАМ 4907 14711 
СОМБОТ, (20) 4900,1 14691 
СОМЪОГ, (30) 4932,1 14782 

















Сравнение результатов показывает, что осесимметричные модели согласуются с очень высокой точностью. 
Разница между трехмерной моделью и осесимметричной может быть вызвана сложностью построения объемной сет- 
ки на наиболее тонких участках геометрии. 

Возможности постпроцессора пакета СОМ$ОГ, позволяют построить трехмерное изображение по итогам ре- 
шения осесимметричной задачи. Здесь и далее в работе приводятся результаты постпроцессорной обработки решений 
осесимметричной задачи. 

Была проведена оценка влияния геометрических параметров на собственные частоты (СЧ). На 
рис. 3 представлены зависимости СЧ от различных геометрических параметров. 
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Рис. 3. Влияние на собственную частоту геометрических параметров: толщины В пьезопреобразователя (а); параметра кривизны К 
направляющей концентратора колебаний (Ъ); толщины 6 концентратора колебаний (с); длины О стержня (4) 


Ею. 3. Ефес! оп авеп/гедиепсу о} зеотейлс рагатаету: цискпезх$ й оррлего@есй“с гапзюгтег (а); сигуаиге рагатеег К оГозсШайоп 
сопсетгейог эшае (Ъ); ИисКпе5з Ь оГозсШайоп сопсетгайот (с); го4 (епай О (4) 
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Из рисунков видно, что с увеличением толщины преобразователя, равно как и с увеличением длины стержне- 
вой части, первая собственная частота уменьшается. Это уменьшение будет значительным, когда приходится увели- 
чивать длину стержневой части, но практически отсутствует при малом изменении толщины преобразователя. Харак- 
терно, что на собственную частоту также практически не влияет кривизна концентратора напряжений и его длина. 

Другой важной характеристикой является амплитудно-частотная характеристика (АЧХ) продольных колеба- 
ний вблизи резонансных частот. Так, на рис. 4 приведены результаты расчета АЧХ продольных колебаний вблизи ре- 
зонансной частоты при определенных значениях независимых физических и геометрических параметров. 
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Рис. 4. АЧХ вынужденных колебаний скальпеля при различных значениях вязкости 
Е. 4. ЕВЕ оГрюгсе4 зсаре! озсШаноп5 ай уатои5 уашез ору15созйу 


Расчеты проводились с учетом акустической жидкости, в которую вставлен конец стержневой части. Воздей- 
ствие акустической жидкости приводит к демпфированию, при котором частота резонанса меняется в зависимости от 
свойств среды. На рис. 4 изображены АЧХ продольных колебаний для различных значений динамической вязкости. 
Второй коэффициент вязкости принимался равным нулю. 

Интересно отметить, что при изменениях вязкости акустической среды в большом диапазоне АЧХ изменяется 
незначительно. Отсюда можно сделать вывод, что основным фактором, влияющим на амплитуду колебаний в окрест- 
ности первого резонанса продольных колебаний, является затухание в упругой среде. Таким образом, увеличение вяз- 
кости акустической среды более чем на порядок не приводит к значительному снижению амплитуды колебаний 
стержневого звена. Интересно также отметить, что вязкость оказывает влияние на амплитуду колебаний, но не на ча- 
стоту. Это может быть связано с небольшой площадью сечения скальпеля. 

Обсуждение и заключения. В настоящей работе были исследованы резонансные частоты продольных 
колебаний системы, состоящей из цепочки разных по физическим параметрам элементов: пьезокерамического, двух 
упругих и акустического. С этой целью были построены осесимметричные и трехмерные конечноэлементные модели. 
Для сравнения были выбраны два конечноэлементных пакета СОМ$ОГ, и АСЕГАМ. Обнаружено хорошее совпадение 
результатов по собственным частотам для осесимметричной задачи. Построенная трехмерная модель послужит в 
дальнейшем при моделировании несимметричных задач. Проведены оценки влияния геометрических параметров на 
значение первой собственной частоты. В частности, рассмотрены варианты, в которых варьировались следующие па- 
раметры: толщина пьезопреобразователя; кривизна и толщина концентратора колебаний; длина стержневого звена. 
Другим важным параметром исследования явилась АЧХ продольных колебаний вблизи первой резонансной частоты. 
Отмечено, что малое изменение динамической вязкости акустической среды мало влияет на эту характеристику. А 
заметное увеличение этой вязкости приводит лишь к снижению амплитуды колебаний, но практически не изменяет 
частоту резонанса. Полученные результаты могут представлять интерес при конструировании ультразвуковых режу- 
щих медицинских приборов. 
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